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ABSTRACT 


A numerical  simulation  of  the  two- 
dimensional  shear  layer  occurring  at  the  interface 
of  two  streams  of  different  velocity  is 
described.  The  shear  layer  is  represented 
by  a large  number  of  discrete  vortices  and  is 
allowed  to  evolve  in  accordance  with  the  velocity 
potential  of  the  vortices.  It  is  found  that  the 
model  effectively  simulates  the  large-scale 
evolution  of  the  layer.  The  model,  however,  is 
found  to  be  inadequate  for  predicting  the  hydro- 
dynamic  sound  generated  in  the  shear  layer 
because  small-scale  vortex  motion  plays  a major 
role  in  acoustic  prediction. 


INTRODUCTION 

BACKGROUND 

The  two-dimensional  mixing  layer  generated  at  the  interface  of  two 

streams  having  different  velocity  and  not  necessarily  having  the  same 

physical  properties  has  recently  been  the  subject  of  much  theoretical  and 
1* 

experimental  study.  This  layer  is  an  archetype  of  all  shear  flows 
because  it  demonstrates  in  simple  fashion  the  processes  of  turbulent  energy 
production,  mixing,  and  entrainment  which  are  fundamental  to  such  flows. 
Much  of  the  research  to  date  has  been  devoted  to  examination  of  the  "large 
coherent  structure"  of  the  turbulence  within  the  layer.  It  appears  that 
the  large  scale  turbulent  motion  within  the  layer,  and  in  other  turbulent 
shear  flows,  is  in  some  sense  deterministic.  The  individual  events 
associated  with  this  large  scale  motion  (for  example,  eddy  growth  and 
interaction)  occur  in  a particular  sequence  and  reoccur  randomly,  but  with 
a statistically  definable  mean  period.  Of  course  the  importance  of  large 


*A  complete  listing  of  references  is  given  on  page  31. 
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turbulent  eddies  was  postulated  by  Townsend  over  twenty  years  ago,  but 

these  eddies  were  viewed  as  entirely  random,  short-lived  phenomena.  However, 

the  visual  and  quantitative  observations  of  the  development  of  spatially 

3 

coherent  patterns  within  the  turbulent  mixing  layer  by  Brown  and  Roshko 

4 

and  by  Winant  and  Browand,  and  within  a turbulent  boundary  layer  by  the 

5 6 7 

Stanford  University  group,  ’ ’ have  brought  about  a redirection  of  much 
research  into  the  turbulence  problem  and  a revision  of  many  long  accepted 
views  as  to  its  nature.  The  current  view  is  that  the  large  scale  motion, 
in  particular  the  interaction  between  the  coherent  structures,  dominates 
the  evolution  and  growth  of  the  shear  layer.  Similar  importance  has  been 
attached  to  large  structures  observed  in  boundary-layer  flows  and  in  the 
near-nozzle  mixing  region  of  axisymmetric  jets.®’^’ 

The  new  awareness  of  the  significance  of  large  scale  coherent 
structures  in  turbulent  flows  has  led  several  researchers  to  suggest  that 
these  structures  might  play  a role  in  the  not  well  understood  mechanisms 
responsible  for  the  hydrodynamic  sound  generated  in  shear  flows. ^ If  this 
supposition  is  correct  and  if  such  structures  are  in  some  sense  determin- 
istic, then  major  advances  in  the  prediction  of  hydrodynamic  sound  might  be 
within  reach. 

The  significance  and  role  of  large  structures  in  sound  generation,  is, 

however,  a subject  of  much  controversy.  For  example,  the  experiments  of 
8 9 10 

Crow  and  Champagne,  Moore,  and  Lau  and  Fisher  have  amply  demonstrated 
that  large  structures  are  a major  if  not  dominant  flow  constituent  in  the 
first  four  to  six  diameters  downstream  of  the  exit  of  an  axisymmetric  jet, 
and  it  is  a widely  held  belief  that  at  least  half  of  the  total  sound  power 
emitted  by  such  jets  originates  from  sources  within  the  first  four 
diameters.  A connection  between  the  coherent  structures  and  at  least  part 
of  the  acoustic  output  might  thus  be  expected.  Crow's  experiments,^  in 
fact,  suggest  a close  relationship.  Crow  performed  far-field  acoustic 
experiments  in  which  a jet  was  excited  internally  over  a range  of  Strouhal 
numbers  (St  = fD/Ug  where  f is  the  excitation  frequency,  D is  the  nozzle 
diameter,  and  Ug  is  the  exit  velocity).  He  measured  a peak  gain  of  as 
much  as  30  dB  (or  a factor  of  about  2500)  in  the  far-field  noise  filtered 
about  a Strouhal  number  of  0.3.  This  Strouhal  number  corresponded  to  a 
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preferred  frequency  for  the  growth  of  disturbances  in  the  jet  mixing  layer 

and  for  the  generation  of  eddy-like  large  structures.  Crow  found  no  change 

in  the  broad-band  noise  in  response  this  forcing.  In  i more  recent  but 

9 

quite  similar  study,  however,  Moore  found  that  even  with  excitation  at  a 

preferred  frequency  for  disturbance  growth  (St  % 0.5  in  his  experiment) 

there  was  no  significant  increase  in  the  far-field  noise  filtered  about  this 

frequency.  There  was,  moreover,  an  increase  of  5-8  dB  in  broadband  noise. 

The  results  of  Crow  and  Moore  are  strongly  at  variance  and,  as  Crighton 
12 

notes,  even  the  fact  that  Crow's  exit  velocities  were  slightly  supersonic 
and  Moore's  subsonic  does  not  explain  the  differences.  Crow's  results 
suggest  that  the  large  eddy  behavior  is  directly  responsible  for  significant 
acoustic  output,  whereas  Moore  has  suggested  that  the  sound  is  generated  by 
smaller  scale  turbulent  motion  which  may  be  modified  through  a highly  non- 
linear process  by  the  excitation  of  large  scale  eddies. 

Although  the  experimental  results  of  Crow  and  Moore  differ 
substantially,  both  indicate  that  large  coherent  structures  play  a signif- 
icant noise-producing  role.  There  have  been  several  attempts  to  numerically 
predict  the  hydrodynamic  sound  generated  in  shear  flows.  In  all  of  these 

the  assumption  that  large  scale  structure  plays  a dominant  role  was 

1 3 

implicit.  Davies  et  al.  attempted  to  develop  a model  based  on  a discrete 

vortex  representation  of  the  mixing  layer  for  the  calculation  of  the  noise 

from  an  axisymmetric  jet.  Although  a peak  in  the  computed  noise  spectrum 

was  found  at  about  St  = 0.5,  most  of  the  characteristics  of  jet  noise  were 

1 4 

not  well  represented.  Hardin  and  Mason  developed  a discrete  vortex 

model  for  noise  generated  by  shear  flow  over  the  mouth  of  a cavity.  While 

they  include  in  their  paper  no  direct  comparison  with  experiment,  they  do 

claim  that  the  model  explains  certain  frequency  and  directivity  phenomena. 

Neither  Davies  et  al.  nor  Hardin  and  Mason  included  mean-flow  refraction  or 

"shrouding"  effects.  In  unpublished  work  (1975),  Lugt  et  al.  attempted  to 

calculate  the  noise  generated  in  a two-dimensional  mixing  layer  represented 

by  discrete  vortices.  Using  no  more  than  48  discrete  vortices,  they  found 

that  they  could  not  adequately  represent  the  mixing  layer  and  hence  that 

the  calculated  sound  power  had  no  physical  meaning.  Both  Hardin  and  Mason 

1 5 

and  Lugt  et  al.  based  their  models  for  noise  generation  on  Powell's 
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formulation  of  the  acoustic  source  terms.  Using  a more  direct  approach, 
Metcalfe  and  Orszag1^  attempted  to  model  jet  fluctuations  with  solutions 
of  the  complete  three-dimensional  Navier-Stokes  equations.  According  to 
Metcalfe^  this  effort  was  not  successful  due  mostly  to  an  inability  to 
obtain  sufficient  statistical  properties  from  the  flow  simulation. 

This  report  presents  the  results  of  an  effort  to  develop  a mathematical/ 
numerical  model  for  the  acoustic  radiation  of  a simple  two-dimensional 
mixing  layer  subject  to  an  applied  disturbance.  It  examines  discrete 
vortex-type  approximations  for  the  prediction  of  hydrodynamic  sound  and 
hopefully  provides  some  insight  into  their  utility. 

MIXING  LAYER  SIMULATION  - LARGE  COHERENT  STRUCTURES 

Several  recent  attempts  to  model  the  behavior  of  the  mixing  layer 

have  emphasized  the  characteristics  of  the  large  scale  structure.  Winant 
4 

and  Browand  proposed  a very  simple  model  which  consists  of  two  vertically 
offset  infinite  arrays  of  equispaced  line  vortices.  According  to  the 
authors  the  vortices  were  found  to  rotate  about  one  another  and  "draw 
closer  together"  (presumably  in  the  mean-flow  direction)  while  the  distance 
between  the  vortex  rows  increased.  Qualitative  agreement  between 
predictions  of  the  model  and  experiment  was  asserted  for  the  layer  growth 
rate. 

1 8 

Acton,  however,  argued  that  a model  for  the  evolution  of  the  mixing 
layer,  especially  for  the  coalescense  of  large  vortical  structures,  requires 
consideration  of  finer  details  of  the  layer  structure.  Acton  therefore 
modeled  the  mixing  layer  with  a distribution  of  discrete  vortices  such  that 
a single  coherent  structure  consisted  of  an  amalgam  of  many  discrete 
vortices.  He  considered  a spatially  periodic  disturbance  of  an  infinite 
discretized  shear  layer  such  that  each  period  consisted  of  two  sinusoidal 
"waves"  displaced  vertically  relative  to  each  other.  Initially,  vortices 
in  each  of  the  "waves"  rol  1 up  into  rotating  concentrations  which  subse- 
quently revolve  about  each  other  until  they  coalesce  (see  Figure  2 of  this 
report  for  the  results  of  a similar  simulation).  The  behavior  reported  by 
Acton  is  in  good  qualitative  agreement  with  the  "pairing"  description  of 
Winant  and  Browand. 
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Most  recently,  Ashurst  reported  a numerical  simulation  of  the 

mixing  layer  using  the  discrete  vortex  approach  without  the  assumption  of 

periodicity.  He  was  able  to  obtain  qood  agreement  with  the  experimental 

20 

data  of  Browand  and  Weidman  for  the  mixing  layer  growth  and  the  time- 
averaged  Reynolds  stress  profile  within  the  layer.  These  calculations  used 
several  thousand  discrete  vortices,  were  quite  expensive,  and  yielded  a 
very  fine  color  film  of  the  shear-layer  evolution. 

The  research  reported  here  extends  the  basic  approach  of  Acton, 
incorporates  additional  vortices,  and  attempts  to  compute  from  the  vortex 
motion  the  hydrodynamic  sound  generated  by  the  shear-layer.  The  assumption 
of  periodicity  is  not  overly  restrictive  for  the  purpose  of  this  study,  and 
it  allows  the  major  aspects  of  shear  layer  behavior  to  be  simulated 
economical ly. 

DISCRETE  VORTEX  SHEAR-LAYER  SIMULATION 

We  consider  a two-dimensional  shear  layer  which  extends  infinitely 
far  in  the  positive  and  negative  x-j -directions  (Figure  1).  The  flow 
velocity  far  above  the  layer  is  -U  and  the  velocity  below  is  U.  The  shear 
layer  of  thickness  6g  at  time  tg  represents  the  lateral  region  (in  x^)  over 
which  the  velocity  varies  from  -U  to  U.  The  velocity  components  in  the 
x-|-  and  X2-directions  are  denoted  by  u-j  and  u^,  respectively.  In  order  to 
investigate  the  nonlinear  evolution  of  the  shear  layer  which  results  from 
the  imposition  of  a spatial  disturbance,  we  represent  the  vorticity 
u>  = sup/sx-i  - 9 u i / 9X2 , within  the  layer,  by  a distribution  of  discrete 
point  vortices  each  with  circulation  r.  We  distribute  the  vortices  in 
cyclic  fashion  over  longitudinal  sections  of  length  A with  N vortices  per 
section,  and  we  consider  the  motion  of  the  vortices  within  one  cycle  and 
assume  that  the  motion  in  all  other  cycles  is  identical.  The  motion  is 
thus  periodic  with  period  A.  We  assume  that  the  effects  of  viscosity  are 
important  only  at  very  small  scales  of  motion  so  that  they  may  be  safely 
ignored. 


EQUATIONS  OF  MOTION 

The  complex  potential  w = u-j  + iu2  of  a single  infinite  row  of  equi- 
distant vortices  of  strength  r and  spacing  x is 

"<2>  ■ S s1"(t) 

where  z = X|  + i x^.  Separating  the  real  and  imaginary  components  of  w,  we 
obtain 

r sinh  2 it  x^/A 

u1(x],x2)  - 2x  cosh  2tt  x2/a  - COS  2tt  x^X 

and 

r sin  2tt  x^/x 

u2(x1’x2)  = 2x  cosh  2tt  x2/X  - cos  Zn  x^/x 

Each  of  the  N vortices  distributed  within  a cycle  may  be  viewed  as  a 
member  of  a different  infinite  row  so  that  we  represent  the  layer  by  N rows 
of 'equidistant  vortices.  Then  the  velocity  at  any  point  induced  by  al  ; N 
rows  of  vortices  is  given  by 

N.  / r ^ sinh  2h(x2-x2  .)/X 

ul(xl’x2)  “ ^ ~\Z\)  cosh  2tt(x2-x2  j)/x-cos  2tt(x1-x1  j)/x  ^ 

and 

N_  f p ^ sin  2ir(x1-x1  .)/X 

u2(xi  ,x2)  = X (jx)  cosh  2tt(x9-x’J  J/X-COS  2n(x,-X,  j/x  ^ 

where  (x,  . ,x„  .)  is  the  position  of  the  jth  vortex  within  the  cycle.  Note 
I » J ^ » J 

that  as  x2  approaches  ±°°,  u^  approaches  + U = rN/2x  so  that  the  required 
vortex  strength  r is  2xU/N. 

In  inviscid  incompressible  flow  vorticity  convects  with  fluid 
elements  so  that  Equations  (1)  and  (2)  can  be  used  to  describe  the  motion 
of  each  vortex,  i .e. , 


d x 


X1  i 

L»  1-  = 11  t 

dt  V 


xl,j,x2,j)  E Ul,j 


(3) 


and 
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I 


d x 


1^"  " U2(xl  ,j,x2,j)  5 U2 , j 


(4) 


Equations  (3)  and  (4)  determine  the  motion  of  the  jth  vortex  in  terms  of 
the  N-l  other  vortices  in  the  cycle.  The  velocities  u,  . and  u9  . are 
given  by  Equations  (1)  and  (2).  For  purposes  of  computation,  we  non- 
dimensional  ize  as  follows: 

uf  = 2xu]/rN 
u*  = 2Xu2/rN 
X*  “ X-|  /x 
x|  = x2/x 

and  t*  = trN/2x2 

Then 


dx^ 

~dt* 


LJ-  = 


1 >i 


1 N sinh  2n(x*  .. -x|  ) 

" j?i  cosh  2'>(xfi'1-xf;j)-cos  2i>(xf_rx^j) 
j7i 


(5) 


and 


dxl,i 

dt* 


1 


N sin  2tt ( x?  . -x*  .) 
y 1 ,i  1 ,3 


2,i  ' N jf-i  cosh  2n(x*ji-x|;j)-cos  2,(x*J-x*J) 


(6) 


and  u^  approaches  +1  as  x*  approaches  + °°  . 


Equations  (5)  and  (6)  constitute  a system  of  2N  coupled,  non-linear, 
ordinary  differential  equations.  With  suitable  initial  conditions  (the 
x.j  ,x2  position  of  each  vortex)  these  equations  determine  the  subsequent 
motion  of  the  vortices  and  the  evolution  of  the  layer.  In  this  study  the 
equations  were  solved  numerically  using,  during  the  early  stages,  a fourth- 
order  Runge-Kutta  method;  solutions  from  which  acoustic  radiation  were 
calculated  were  obtained  with  a sixth-order  Milne  predictor-corrector 
algorithm  since  it  was  found  that  much  greater  accuracy  was  required. 
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All  the  calculations  described  in  this  report  were  performed  on  the 
Texas  Instruments  Advanced  Scientific  Computer  at  the  Naval  Research 
Laboratory  with  double  precision  arithmetic.  DTNSRDC  users  access  this 
computer  from  a remote  job  entry  terminal  at  DTNSRDC.  Certain  data  were 
transferred  via  punched  cards  and  magnetic  tape  to  the  CDC  6700  at  DTNSRDC 
where  the  TIMAGE  interactive  graphics  system  and  the  DISSPLA  computer 
graphics  package  were  used  to  generate  most  of  the  figures. 

SHEAR-LAYER  EVOLUTION 

Figures  2a  through  2h  show  a sample  simulation  of  the  evolution  of 
a shear  layer  represented  by  256  vortices  per  period  x.  The  shear  layer  of 
initial  thickness  0.05  is  assumed  to  have  been  disturbed  by  some  external 
means  and  at  t=0  has  the  configuration  shown  in  Figure  2a.  Note  that  each 
wavelength  consists  of  two  sinusoidal  "waves",  each  of  length  1/2,  which  are 
displaced  vertically  relative  to  each  other  a distance  0.02.  Each  of  the 
waves  has  amplitude  -0.05.  Vortices  in  the  left-hand  and  right-hand  "waves" 
are  marked  with  0's  and  X's,  respectively.  The  integration  was  performed 
with  the  Runge-Kutta  method  with  time  step  0.005  and  required  101  CP  seconds 
to  reach  t*=1.0.  During  the  initial  stages  of  the  evolution  (Figures  2a-2d), 
the  vorticity  in  both  "waves"  concentrates  in  nearly  identical  fashion. 

(The  behavior  would  be  identical  were  it  not  for  the  vertical  displacement.) 
As  the  shear  layer  continues  to  evolve,  the  two  vortex  amalgams  interact, 
rotate  about  each  other  (Figures  2e-2g),  and  coalesce  (Figure  2h).  This 
process  is  quite  similar  to  the  "pairing"  process  observed  experimentally 

4 

by  Winant  and  Browand.  The  elongation  of  the  structures  in  the  flow  (x^) 
direction  (see  especially  Figure  2e)  was  also  observed  experimentally  by 
Winant  and  Browand  who  described  the  vortices  as  taking  on  a "more  or  less 
elliptic  shape  with  major-to-minor  axis  ratio  of  2."  It  is  also  clear  that 
fluid  is  transported  to  the  opposite  side  of  the  mixing  layer  before  being 
ingested  into  the  central  region.  Again  similar  behavior  was  noted  by 
Winant  and  Browand. 

On  the  basis  of  this  simulation  and  other  simulations  we  have  per- 
formed, we  believe  that  discrete  vortex  simulation  can  effectively  model  the 
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large  scale  behavior  of  the  shear  layer.  Acton,  who  used  only  96 
vortices,  performed  extensive  calculations  with  his  model  and  has  presented 
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Figure  2 - Evolution  of  a Shear  Layer  Represented  by  Discrete 
Vortices  Shown  at  Different  Times 
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Figure  2 (Continued) 
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his  results  in  great  detail,  including  comparisons  with  experimental  data. 
The  overall  characteristics  of  the  shear-layer  evolution  shown  in  Figures 
10a  through  lOh  in  his  paper  are  in  excellent  agreement  with  the  calculation 
presented  above.  His  paper  is  recommended  to  the  interested  reader. 


VORTICITY-GENERATED  HYDRODYNAMIC  SOUND 
MATHEMATICAL  FORMULATION 

To  develop  an  estimate  for  the  acoustic  radiation  from  the  evolving 

21 

shear  layer  we  begin  with  Lighthill's  fundamental  expression  for  the 
acoustic  radiation  generated  by  a compact  region  of  fluctuating  flow  in  an 
otherwise  quiescent  fluid: 


^ ' C0  ■’■’■I  <7> 

where  I % Pg'/  ; v , c^  and  are  the  ambient  undisturbed  sound  speed  and 
density,  and  v is  the  fluid  velocity  computed  with  the  assumption  of 
incompressibility.  In  applying  Equation  (7)  to  the  shear  layer  we  assume 
that  U/Cq  <<  1.  Powell  showed  that  the  origin  of  the  sound  field  can  be 

associated  with  the  fluctuating  vorticity  field;  another  derivation  was 

22  2T 

given  by  Howe.  We  note 

VV-T  = pQ  v-L  + p0v2(1v2  ) (8) 

-y 

where  L is  defined  by  L = u x v and  w is  the  local  vorticity.  Thus, 
specification  of  L = u x v over  the  fluctuating  region  (which  is  assumed 
to  be  small--its  characteristic  dimension  is  much  smaller  than  an  acoustic 
wavelength)  yields  an  estimate  for  the  acoustic  radiation. 

The  formal  solution  to  Equation  (7)  with  the  source  term  of 
Expression  (8)  is  obtained  using  a Green's  function  and  is 
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p(3t»t)-pn  , r r + 6(t-T-|x-?|/r) 

= 4^  //(VL>  “ d3J  d’ 


|x-y| 


kllHr) 


4itcq 


2\  6(t-T-|x-y|/c0) 

— d3y  dx 


|x-y| 


This  kind  of  solution  is  sometimes  termed  a retarded  potential.  The 

acoustic  signal  reaching  x at  time  t originates  at  position  y within  the 

fluctuating  region  at  t-|x-y|/cn.  We  apply  the  usual  simplifications  to 
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this  integral  (for  example,  see  Lighthill  and  Crighton  ).  We  make  use 
of  the  properties  of  linear  differential  operators  applied  to  convolution 
products  and  then  perform  the  x-integration  to  obtain 


— rr'x  ' / L (y.t-|Sf-y|/c  ) 

P0  4"cn  i 0 |x-y| 


* Jr 


d3y 


|x -y| 


We  then" perform  the  x-differentiations  and  ignore  higher -order  terms 
so  that 


p(x,t)-p 


0 . 


4ttC?  I ( 


. ) -^(y.t-|3-y|/c0)d3> 

o ? c0>x-yi2 


+ — / 
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32  V2  ,+ 


4«5  j cq| x - y | at2 


\ (y,t-|x-y|/c0)d3y 


Finally,  we  assume  that  the  eddy  region  is  compact,  that  is,  its 
characteristic  linear  dimension  is  much  smaller  than  the  acoustic  wave- 
lengths of  the  sound  it  generates,  and  that  x-y  % x,  so  that  we  consider 
only  the  far-field  acoustic  radiation.  We  obtain 
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L(y,t-|x|/c0)d3y 
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The  last  term  may  be  neglected  for  compact  sources  so  that 


p(x,t)-p. 


4*cj*|x|  3t2  y 


TTT  • L(y,t-|x|/cQ)d3y 
x u 


Equation  (9)  shows  that,  to  first  order,  the  acoustic  radiation  is 

•f 

determined  from  the  second  time  derivative  of  L.  Note  that  the  acoustic 
signal  at  x at  time  t is  now  related,  through  the  second  time  derivative, 
to  fluctuations  at  the  same  time,  t-|x|/cn,  for  all  locations  y. 

We  may  evaluate  L and  its  time  derivatives  in  a simple  fashion  for  an 
array  of  discrete  point  vortices.  For  such  an  array 


w(y.t)  - I 
j=l 


r 6(xl,j'yl^6(x2,j"y2^e3 


where  e3  is  a unit  vector  normal  to  the  plane  of  motion.  We  substitute 
Equation  (10)  into  Equation  (9)  and  perform  the  y integration  (which  is 
made  easy  by  the  appearance  of  6 functions)  to  yield  the  acoustic  radiation 
per  unit  length  (in  the  e3-direction)  at  x=x^e1+x2e2+x3e3 


P(x,t)-p 


o = L 

4ttcI?  I 


dt2  i*! 


j.i  \ i*i  / 


(re3  x5j) 


where  1).  = x . , the  velocity  of  the  jth  vector,  and  u.  and  x.  are  evaluatec 
J J J J 

at  t- | x | /Cq . This  acoustic  radiation  is  due  to  flow  fluctuations  in  the 

plane  x3  = constant.  Equation  (11)  is  written  in  terms  of  the  previously 

defined  non-dimensional  variables  as 
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p0A^3 
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2it|x*|N  dt*2  j = 


N /x*  • x*.\ 
j?iV  | x*  | / 


(e3xu*) 


(12) 


where  M = tN/2aCq  = U/Cq  is  the  mean  flow  Mach  number,  and  Ay*  = Ay3/A. 
We  may  express  Equation  (12)  in  terms  of  polar  coordinates  r*  = |x*|  and 
0 = tan~^x£/x|  in  the  plane  x^  = constant  as 

p ( r* , 0 , t*)  ~Pq  Iy|4 


p0ay3 


2nr* 


[(A-B)cos  20  + C sin  20  + A + B] 


(13) 


where 


= _ 1 v 
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and 


B = 2N  p x2 , j u*,j 

l N 

r — — ( x*  ij*  - x*  u*  ^ 

2N  jt,  (l,j  Ul,j  X2,j  U2 , j 


The  quantities  A,  B,  and  C are  evaluated  at  time  t*-M|x*|. 

CALCULATION  OF  SHEAR-LAYER  HYDRODYNAMIC  SOUND 

Equation  (13)  provides  a means  for  estimating  the  hydrodynamic 
sound  generated  by  an  unsteady  vortical  flow  region  approximated  by  a 
distribution  of  discrete  point  vortices.  In  this  section  we  describe  our 
application  of  Equation  (13)  to  the  sound  emitted  by  a single  spatial  period 
of  the  shear  layer. 

The  presence  of  third  time  derivatives  of  position  in  the  expression 
poses  significant  difficulties.  We  found  that  numerical  estimation  of 
these  derivatives  with  finite  differences  did  not  produce  sufficiently 
accurate  results  without  severe  restrictions  on  the  time-step  size.  We 
therefore  chose  to  obtain  these  derivatives  by  analytically  differentiating 
the  equations  of  vortex  motion.  Equations  (5)  and  (6).  The  resulting 
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expressions  directly  relate  x,  .,  x9  .,  x,  .,  and  iu  . to  x,  . and  x9  . 

without  truncation  error  (except  that  introduced  in  the  calculation  of  x,  . 

* »J 

and  x,  ■).  Note  that  all  of  the  shear  layer  simulation  results  reported 
in  this  section  were  obtained  with  a sixth-order  Milne  predictor-corrector 
integration  scheme  applied  to  Equations  (5)  and  (6). 

Figures  3a  and  3b  display  two  stages  in  the  evolution  of  the  vortex 
system  in  the  test  case  of  motion  with  no  external  shear-layer  disturbance. 
The  layer  is  approximated  by  126  vortices  per  unit  length  distributed 
initially  in  18  equispaced  columns,  7 vortices  per  column,  and  the  calcu- 
lation was  performed  with  At*=0. 00125.  This  smaller  number  of  vortices  was 
used  in  acoustic  model  evaluation  calculations  to  reduce  the  otherwise  very 
large  expense  of  performing  computations  with  very  small  time  steps.  The 
figures  suggest  a great  deal  of  energetic  non-physical  (numerical)  vortex 
activity,  and  by  t*=0.5  (Figure  3b)  significant  numerical  round-off  error 
has  accumulated.  The  regular  pattern  along  the  layer  and  the  antisymmetry 
across  the  layer  of  vortex  positions,  which  we  expect  in  an  exact 
calculation,  is  significantly  eroded.  We  emphasize  that  this  error  is  due 
to  accumulated  round-off  in  the  calculation  of  the  N-term  sums  in 
Equations  (5)  and  (6),  and  not  to  truncation  error  in  the  integration 

l 

scheme.  Truncation  error  affects  the  position  of  the  vortices  in  each  row 
in  an  identical  manner  so  that  a regular  vortex  pattern  would  still  result. 

The  acoustic  radiation  as  predicted  from  Equation  (13)  for  the 
simulation  described  here  is  displayed  in  Figures  4a  and  4b  as  a function 
of  angle  9 measured  from  the  position  x^-direction.  The  radiation  is 
presented  in  terms  of  a normalized  "instantaneous"  intensity  which  we 
define  as 


..  ..  ..  ..  ..  2 

and  which  is  equal  by  Equation  (13)  to  [(A-B)cos  20  + C sin  2e  + A + B]  . 

Dividing  this  normalized  intensity  computed  at  time  t*  by  (2irr*)2  and 

multiplying  by  M8  yields  [(p-Pg)/(pgAy*)]2 , which  is  proportional  to  the 

"instantaneous"  intensity,  a^p2/pQ,  of  the  sound  at  position  r*,  e at  time 
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Figure  3 - Position  of  Vortices  in  Shear  Layer 
With  No  Applied  Disturbance 
(At  t*=0.0,  vortices  are  arranged  in 
columns,  seven  per  column) 
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t*  + Mr*.+  Also  shown  in  Figures  4 is  the  intensity  predicted  by  a second 
identical  simulation  with  At*=0.0025.  At  t*=0.25  and  at  earlier  times  the 
intensities  predicted  by  the  two  calculations  are  in  excellent  agreement,  and 
both  the  vortex  positions  and  the  instantaneous  vortex  velocities  were  found 
to  agree  to  within  more  than  four  significant  figures.  At  t*=0.5,  however, 
significant  differences  in  the  predictions  became  apparent.  These  differ- 
ences are  caused  by  very  slight  differences  in  vortex  position  between  the 
two  calculations  which  result  both  from  the  round-off  error  already  mentioned 
and  truncation  errors.  These  differences  in  position  lead  to  large  differ- 
ences in  the  vortex  velocities,  accelerations,  and,  hence  the  acoustic  out- 
put. Since  the  discrete  vortex  approximation  is  made  with  the  assumption 
that  the  fine  details  of  the  flow  are  not  of  consequence,  this  result  throws 
doubt  on  the  applicability  of  the  model  for  predicting  hydrodynamic  noise. 

The  significant  dependence  of  the  acoustic  radiation  on  the  detailed 
vortex  distribution  at  each  instant  appears  in  more  dramatic  fashion  in  a 
second  series  of  calculations.  Figures  5a  through  5c  and  6a  through  6c 
display  vortex  positions  at  three  time  levels  from  two  calculations  of  the 
layer  evolution  resulting  from  a single-wave  disturbance.  In  the  simula- 
tion depicted  in  Figures  5a  through  5c  126  vortices  are  "stacked"  in 
18  columns,  7 vortices  per  column.  In  the  second  simulation  (Figures  6a 
through  6c)  126  vortices  are  "staggered"  in  18  columns  with  4 vortices  per 
column,  and  18  columns  with  3 vortices  per  column.  Although  the  over-all 
evolution  of  the  shear  layer  as  represented  by  the  two  computations  is 
quite  similar,  the  details  of  the  vortex  distributions  are  very  different 
as  might  be  expected.  The  differences  in  the  details  of  the  vortex 
distributions  have,  unfortunately,  a very  significant  effe:t  on  the  acoustic 
radiation.  (Figures  7a  through  7c.)  The  predicted  intensities  are  quite 
different  both  in  magnitude  and  directivity.  These  results  indicate  that 
the  hydrodynamic  sound  computed  from  the  discrete  vortex  model  is  not 
dominated  or  controlled  by  the  large  scale  evolution  of  the  vortex  system, 
but  is  in  fact  strongly  dominated  by  high  frequency  components  which  are 

The  intensity  of  sound,  I,  is  defined  as  the  average  rate  at  which  energy 
is  transported  across  a unit  area  normal  to  the  propagation  direction,  so 
that  I = Cq<p2>/pq  where  < > denotes  time  average. 
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Figure  5c  - t*=0.5 

Figure  5 - Position  of  Vortices  in  Shear  Layer  Resulting 
From  a Single-Wave  Disturbance 
(At  t*=0.0,  vortices  are  arranged  in  columns, 
seven  per  column) 
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Figure  6 - Position  of  Vortices  in  Shear  Layer  Resulting 
From  a Single-Wave  Disturbance 
(At  t*=0.0,  vortices  are  arranged  in  columns, 

4 vortices  per  column  and  3 vortices  per  column) 
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Figure  7 - Comparison  of  Normalized  Intensities  of  Acoustic  Radiation  from 
Shear  Layers  at  Three  Different  Times  Resulting  from  a Single-Wave 
Disturbance  and  with  Different  Initial  Distributions  of  Vortices 
(Dashed  lines  correspond  to  the  case  of  the  stacked  initial 
distribution  shown  in  Figure  5a;  the  chain-dotted  lines, 
to  the  staggered  initial  distribution  of  Figure  6a) 


dependent  on  the  detailed  vortex  distribution.  Since  we  cannot  expect  flow 
details  (on  the  order  of  vortex  spacings)  to  be  properly  simulated  with 
discrete  vortices,  the  model  just  described  is  not  appropriate  for  hydro- 
dynamic  sound  calculations. 

The  above  conclusion  is  supported  by  the  earlier  work  of  Lugt  et  al . 
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(1975)  who  used  Hama  and  Burke's  single  row  model  with  12,  24  and  48 
vortices  per  wave  length.  Lugt  found  that  the  acoustic  output  generated  by 
the  single  row  model  was  highly  dependent  on  the  distribution  of  vortices 
which  comprise  the  developing  and  evolving  vorticity  concentration 
(Figure  8).  In  such  simple  models,  however,  it  is  difficult  to  distinguish 
the  effects  of  large  scale  evolution  from  the  effects  of  small  scale  vortex 
interaction. 


OF  SOUND  POWER 


Figure  8 - The  Point  Vortices  Closest  to  the  Center  of  a Rolling-up 
Discontinuity  Sheet  Generate  Most  of  the  Sound.  The  Shaded 
Area  Represents  the  Upper  Side  of  an  Initially  Horizontal 
Discontinuity  Sheet. 

(From  H.J.  Lugt) 

INCLUSION  OF  A FINITE-RADIUS  VORTEX  CORE 

In  order  to  overcome  a number  of  problems  which  have  been  encountered 
in  various  discrete  vortex  flow  simulations,  several  authors  have  intro- 
duced the  concept  of  a finite-radius  vortex  core.  The  assumption  is  made 
that,  at  distances  from  a vortex  less  than  some  specified  value,  the 
velocity  induced  by  that  vortex  is  not  given  by  the  potential  flow  r’1 

distribution  but  is  determined  from  a different  non-singular  distribution. 
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We  incorporated  Chorin's  formulation  of  the  finite-radius  core  into  the 
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discrete-vortex  shear-layer  model.  We  hoped  that  by  smearing  out  the  small- 
scale  details  of  the  layer  evolution  and  thereby  reducing  the  large  acceler- 
ations induced  by  vortices  on  their  close  neighbors,  acoustic  predictions 
less  sensitive  to  the  small-scale  vortex  motion  could  be  obtained. 

Following  Chorin,  we  let  (in  dimensional  form) 

_r_  ^2_ 

U1  " " 2™  r 

and 

_ r X1 

u2  2tio  r 

with  r = /x2+y2,  denote  the  components  of  the  velocity  at  (x1>x2)  induced 
by  a vortex  at  the  origin.  These  expressions  are  applied  at  r < o where  0 
is  the  core  radius;  for  r > 0 the  r'l  profile  is  applied.  The  combined 
velocity  expression  thus  exhibits  the  potential  flow  distribution  outside 
the  core  and  finite  values  of  u-|  and  u^  inside  the  core. 

A number  of  computations  were  performed  with  several  values  of  the 
core  radius  0.  These  calculations  show  that,  even  with  0 as  large  as  one- 
third  of  the  initial  shear-layer  thickness,  the  large-scale  evolution  of  the 
layer  is  not  appreciably  affected.  The  small-scale  behavior  of  the  layer 
is,  however,  drastically  affected  and,  in  general,  the  acoustic  intensity 
is  significantly  lowered.  Figures  9a  through  9c  compare  the  predicted 
intensities  (with  a core  radius  of  one-third  the  initial  shear  layer 
thickness)  for  a computation  beginning  with  the  "stacked"  vortex  distribution 
shown  in  Figure  5a  and  a computation  beginning  with  the  "staggered" 
distribution  shown  in  Figure  6a.  The  intensities  predicted  by  both  the 
"stacked"  and  "staggered"  computations  are  a good  deal  lower  than  those 
predicted  without  the  finite-radius  core  (Figures  7b  and  7c).  This  might 
be  expected  since  the  accelerations  and  time  rates  of  acceleration  change 
of  closely  spaced  vortices  are  substantially  reduced  and  these  quantities 
are  significant  in  the  acoustic  prediction  (see  Equation  (13)).  The 
intensities  at  t*=0.5  predicted  from  the  "stacked"  and  "staggered" 
computations  are  encouragingly  similar  but  by  t*=0.75  they  are  quite 
different  so  that  there  is,  unfortunately,  no  indication  that  the  critical 
influence  of  the  initial  vortex  distributions  on  the  acoustic  radiation 
lessens  as  the  vortex  system  evolves. 
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Figure  9 - Comparison  of  Normalized  Intensities  of  Acoustic  Radiation 
from  Shear  Layers  at  Three  Different  Times  (with  Chorin's  Finite- 
Radius  Core)  Resulting  from  a Single  Wave  Disturbance  and  with 
Different  Initial  Distributions  of  Vortices 
(Dashed  lines  correspond  to  the  case  of  the  stacked  initial 
distribution  shown  in  Figure  5a;  the  chain-dotted  lines, 
to  the  staggered  initial  distribution  of  Figure  6a) 


Figure  9a  - t*=0.25 
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SUMMARY  AND  CONCLUSIONS 


We  have  shown  that  a discrete  vortex  model  of  the  two-dimensional 
shear  layer  effectively  models  the  large  scale  evolution  and  behavior  of  the 
layer  resulting  from  the  application  of  an  imposed  disturbance.  An  identical 

I O 

conclusion  was  reached  by  Acton.  The  over-all  behavior  of  the  layer  is 

in  qualitative  agreement  with  the  experimental  observations  of  Winant  and 
4 

Browand.  However,  although  we  are  aware  of  no  directly  relevant  data 

with  which  to  compare  the  model  estimates  of  the  hydrodynamic  sound 

generated  in  the  layer,  our  investigation  has  led  us  to  conclude  that  the 

vortex  model  is  not  an  appropriate  representation  of  the  shear  layer  for 

acoustic  predictions.  We  find  that  the  detailed  small-scale  motion  of  the 

vortex  system  plays  a major  role  in  the  nonphysical  (numerical)  acoustic 

predictions;  this  result  is  inconsistent  with  the  original  premise  under 

which  all  vortex  models  are  applied;  namely,  that  the  small  scale  behavior 

is  not  of  significance.  Because  small  scale  motion  has  overwhelming 

effects  on  the  acoustic  predictions,  we  are  not  able  to  elucidate  the  role 

of  the  large  scale  behavior  in  the  generation  of  hydrodynamic  noise. 

The  conclusions  presented  here  indicate  that  other  applications  of 
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discrete  vortex  models  to  the  prediction  of  hydrodynamic  noise  ’ should 
be  carefully  scrutinized.  Although  these  other  applications  have  been 
made  to  statistically  steady  flow  fields,  the  time-averaged  results  may 
still  be  significantly  affected  by  the  non-physical  small-scale  vortex 
motion. 
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